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Abstract. This paper deals with exploiting symmetry for solving linear and integer 
programming problems. Basic properties of linear representations of finite groups can 
be used to reduce symmetric linear programming to solving linear programs of lower 
dimension. Combining this approach with knowledge of the geometry of feasible inte- 
ger solutions yields an algorithm for solving highly symmetric integer linear programs 
which only takes time which is linear in the number of constraints and quadratic in the 
dimension. 



1. Introduction 

It is a known fact that many standard (integer) linear programming formulations of relevant 
problems in optimization show a lot of symmetry. In this situation a standard branch-and- 
cut framework repeatedly enumerates symmetric solutions, and sometimes this renders such 
methods useless. To address these issues the last decade saw a number of approaches to 
devise algorithms specialized to symmetric optimization problems. We mention a few: Mar- 
got suggests to solve symmetric integer linear programs (ILPs) via a pruned branch-and-cut 
approach involving techniques from computational group theory [141 115) . A recent improve- 
ment in this direction is "orbital branching" devised by Ostrowski et al. |17j . Friedman [5] as 
well as Kaibel and Pfetsch [5] treat symmetric ILPs by shrinking the domain of feasibility 
by cutting off symmetric solutions. Gatermann and Parrilo apply results from represen- 
tation theory and invariant theory to semidcfinite programming [5], which includes linear 
programming as a special case. Our approach is close in spirit to this paper. See also the 
survey of Margot [13] for a general overview of symmetric integer linear programming. 

This is how our paper is organized: first we analyze linear programs with an arbitrary 
finite group of linear automorphisms. Most results in this section are known. A first key 
observation, Theorem I15[ is that symmetric linear programming can be reduced to linear 
programming over the fixed space of the automorphism group. Sections |3] and [3] translate 
these results to the context of integer linear programming. In the sequel we concentrate 
on groups acting as signed permutations on the standard basis of M". Section [5] contains 
our main contribution: our Core Point Algorithm FBI can solve an integer linear program in 
K" whose group of linear automorphisms contains the alternating group of degree n (acting 
as signed permutations) in 0{mn^) time, where m is the number of constraints. This is 
in sharp contrast with the known NP-completeness of the general integer linear program- 
ming feasibility problem. While our algorithm only works for ILPs with an exceptionally 
high degree of symmetry we believe that this is a first step towards an entire new class 
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of algorithms dealing with symmetry in ILPs. Suitable generalizations are the subject to 
ongoing research. In Section |6] we discuss algorithms to determine groups of automorphisms 
of integer linear programs. This leads to algorithmic problems which turn out to be graph- 
isomorphism-complete. The final Section [5] contains experimental results. One of the ILP 
classes that we investigated computationally is motivated by work of Pokutta and Stauffer 
on lower bounds for Gomory-Clivatal ranks |19) . Section [7] explains the construction. 

We are indebted to Tobias Achterberg, Leo Liberti, Marc Pfetsch, Sebastian Pokutta, 
and Achill Schiirmann for valuable discussions on the subject of this paper. 

2. Automorphisms of Linear Programs 

The purpose of this section is to introduce the notation and to collect basic facts for future 
reference. The results of this section up to and including Corollary [T3] can be found in the 
literature which is why we skip some of the proofs. 
We consider linear programs LP(yl, 6, c) of the form 

, , max c*x 

^ ' s.t. Ax<b, a; e M" 

where A e R"^", b e R", and c e R" \ {0}. Throughout we will assume that the set 
P{A, b) :— {x £ R" I Ax < b} of feasible points is not empty, and hence it is a convex 
polyhedron, which may be bounded or unbounded. We will also assume that an optimal 
solution exists. This is to say, our linear program LP(^, b, c) is bounded even if the feasible 
region may be unbounded. In this case the set of optimal solutions forms a non-empty face 
of P{A,b). Our final assumption c ^ is not essential for the algorithms below, but it 
allows to simplify the exposition somewhat. 

Each row of the matrix A corresponds to one linear inequality. Suppose that one of these 
rows is the zero vector. Then the corresponding right hand side must be non-negative, since 
otherwise the linear program would be infeasible, and this was explicitly excluded above. 
But then this linear inequality is trivially satisfied. Therefore we will further assume that the 
matrix A does not contain any zero rows. In this case each row defines an affine hyperplane. 
This way LP{A,b,c) gives rise to an arrangement ^{A,b) of m labeled affine hyperplanes 
in R". 

Definition 1. An automorphism of the linear program LP(A, 6, c) is a linear transformation 
in GL„R which induces a permutation of 3{{A, 6), which leaves P{A, b) invariant, and which 
does not change the objective value c*x for any feasible point x E P{A, b). 

The objective function is linear, and hence it follows that an automorphism of LP (A, 6, c) 
does not change the objective value on the linear span \m{P{A,b)) of the feasible points. 
One could also take more general affine transformations into account. In all what comes 
below this would require a number of straightforward changes. We refrain from doing so for 
the sake of a clearer exposition. The following examples show that the three properties to 
be satisfied by a linear automorphism are mutually independent. 

Example 2. For rn ~ n — 1 let A ~ 1, 5 = 0, and c — 1. The feasible region is the 
non-positive ray R<o- Multiplication with any positive real number 7 leaves the feasible 
region and the hyperplane arrangement (consisting of the origin) invariant. If 7 7^ 1 the 
objective function is not preserved. 
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Example 3. For m = n = 2 let 

- - -°o ■^-"^ ^- (-„') - 

Then P{A,h) is the non- negative quadrant in R^. Now 7 = (0-1) leaves the coordinate 
hyperplane arrangement 3-C(A, b) invariant, but it changes the feasible region. For each 
X €M? we have c^x = c^'yx. 

Example 4. For m = 3 and n = 2 let 




The feasible region is the non- negative quadrant in M^; the third inequality is redundant. 
The linear transformation 7 = (ij) leaves the feasible region invariant, and it satisfies 
c^x — c*7X for all x S M. The hyperplane arrangement 3-C{A, b) is changed. 

For more examples see [3]. There it is also shown that each finite (permutation) group 
occurs as the group of automorphisms of a linear program. 

Remark 5. It is always possible to scale the rows of the extended matrix {A\b) such that 
the leftmost non-zero coefficient is ±1. This allows to remove duplicate inequalities from 
the input by sorting. The complexity of sorting the rows by pair- wise comparison is of order 
0{mn log m). This can be neglected in the asymptotic analysis of our algorithms below since 
it is always dominated. This way we can always assume that the hyperplanes in 3{.(A,b), 
that is, the inequalities, and the rows of the extended matrix {A\b) are in a one-to-one 
correspondence. In the rational case it is more natural to scale the inequalities to integer 
coefficients which are coprime. This is what we will usually do. For a more sophisticated 
algorithm to sort out equivalent constraints, see Bixby and Wagner [2]. 

Since we view points in R" as column vectors, a matrix 7 representing a linear transfor- 
mation acts by multiplication on the left. The adjoint action on the row space, and thus 
on the set of linear inequalities, is by multiplication of the inverse transpose 7^* on the 
right. The set of linear transformations permuting the arrangement 3{.{A, b) forms a closed 
subgroup of GL„R, that is, a linear Lie group. Similarly, the set of linear transformations 
leaving the feasible region P{A, b) invariant forms a linear Lie group. It follows that the set 
Aut(LP(a, b, c)) of automorphisms of the linear program LP(A, 6, c) also forms a linear Lie 
group. For basic facts about (linear) Lie groups, see Rossmann 22 . 

Remark 6. Clearly, the value and the set of optimal solutions of a linear program only 
depend on the non-redundant constraints. At the expense of one linear program per con- 
straint one can get rid of the redundant ones. This obviously does not help to reduce the 
complexity of solving the linear program given since the linear program for a redundancy 
check is of the same size. However, for more costly algorithmic problems, like integer pro- 
gramming as is discussed below, this reduction can be useful. In particular, this will be 
the case when the group of automorphisms becomes larger, see Example |4l Notice that the 
notion of "invariance" from [6j Definition 3.1], specialized to linear programming, implies 
that redundant constraints are going to be ignored. 

The following result is a consequence of the combinatorial properties of a convex poly- 
tope P: the faces of P are partially ordered by inclusion, and this partially ordered set 
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forms a lattice. The automorphisms of this lattice, the face lattice of P, are called combi- 
natorial automorphisms. Each linear (or afRne or projective) automorphism of P induces 
a combinatorial automorphism, but, in general, a polytope may have many combinatorial 
automorphisms which are not linearly induced. See Ziegler [24i for the details. 

Lemma 7. If the feasible region P{A,b) is bounded and full-dimensional, then the automor- 
phism group Aut(LP(yl, 6, c)) is finite. Moreover, the objective function c satisfies c^jx = c^x 
for all a; e M" and e Aut(LP(^, 6, c)). 

Proof. Let d be a vertex of the polytope P — -P(A, 6). Since dimP = n there are vertices 
witW2, ■ . . each of which shares an edge with and such that the difference vectors wi — 
v,W2 —V, . . . , Wn — V form a basis of R". This implies that each combinatorial automorphism 
of P is induced by at most one linear automorphism. Hence the group Aut(LP(yl, 6, c)) 
is contained in the group of combinatorial automorphisms of P, which is finite. While 
Definition [T] asks that each feasible point is mapped to a (feasible) point with the same 
objective value, the additional claim deals with all points, feasible or not. However, this 
follows from lin(P(yl, b)) — R" as c is linear. □ 

If the polyhedron P{A, b) is not full-dimensional, then the automorphism group is a direct 
product of the group of automorphisms fixing the linear span of P{A, b) with a full general 
linear group of the orthogonal complement. In the sequel we will therefore restrict our 
attention to the fuU-dimcnsional case. 

Definition 8. Given a subset Y C E" and a group F < GLr acting on Y , the set of fixed 
points of Y with respect to an element 7 G F is defined by 



The set of fixed points Fix^(R") is the (possibly zero-dimensional) eigenspace Eig(7, 1) 
of the linear transformation 7 with respect to the eigenvalue 1. This implies that Fixr(R") 
is a linear subspace for any group F of linear transformations. More generally, Fixr(i^) is 
the intersection of this subspace with the set Y . 

Remark 9. If the linear group F < GL„M is generated by the set G C F, then 



In particular, if G is finite, that is, if the group F is finitely generated, this leads to an 
algorithm to compute (a primal or dual basis of) the fixed space by solving one linear 
system of equations per transformation in the generating set G. 

Remark 10. Let F < Aut(LP(yl, 6, c)) be a group of automorphisms of the linear program 
LP(A,5, c) such that P{A,b) is bounded and full-dimensional. Then, by Lemma [71 the 
set of fixed points Fixr(R") contains the one-dimensional linear subspace spanned by the 
objective vector c. 

For any finite set 5" C M" let 



Fix^(y) {yer |7y = y} . 
Therefore, the set of fixed points of Y with respect to F is given by 



Fixr(y) {y e r |7y = y for aU 7 e F} = P|Fix^(y). 

7er 



Fixr(M") = flFix^(M") = f| Eig(7, 1) . 



7GG 7GG 
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be its barycenter. The two subsequent results are basic observations from representation 
theory, not restricted to apphcations in (hnear) optimization. For a brief proof, for instance, 
see [33] Lemma 3.5]. 

Lemma 11. The map 

R" ^ Fixr(]R") , V K> l3{Tv) 
which sends a point to the barycenter of its T- orbit is a linear projection onto the fixed space. 

Let S C R" be a finite set which is spanning, that is, we require hn(S') = R". Further let 
r be a finite subgroup of GL„M acting on S. Phrased differently, we are considering a linear 
representation of an abstract group F on the vector space R" which induces a permutation 
representation on the set S. In this case F splits S into disjoint orbits Oi,02, ■ ■ ■ , Ok- In 
our applications below, S will usually be the set of vertices of some polytope which linearly 
spans R". 

Lemma 12. For the fixed space o/F we have 

Fixr(R") - lin{/3(Oi),/3(02),...,/3(Ofc)}. 
In particular, dimFixr(R") < k. 

Proof. Since S* = Oi U O2 U • ■ • U Ofe is spanning and since the union of the orbits gives S it 
follows that 

(2) R" = lin(Oi) + lin(02) + --- + lin(Ofc). 

For i G {1,2, ... ,k} the linear subspace lin(Oi) is F-invariant. If we apply the surjective 
linear map v n- /3(Fw) from Lemma [TT] to the set S, we obtain a generating set for Fixr(R"). 
Applying the same map to a single orbit Oi similarly yields a generating set for Fixr (lin(O.i)). 
Now the claim follows from the equation FOj = O^. □ 

Notice that the sum decomposition ([2]) is not necessarily direct. We now apply the results 
obtained so far to a finite group of automorphisms of a linear program. 

Proposition 13. Let F < Aut(LP(^, 6, c)) be finite. If x £ R" is an arbitrary point, the 
barycenter of its T-orbit satisfies c^/3(Tx) = c*x. If moreover, x G P{A, b) is feasible, then 
/3{Tx) is feasible, too. 

Geometrically this means that the points of one orbit are in the same affine hyperplane 
orthogonal to c. 

Proof. As the objective function is constant on the orbit Fa; it follows that c*/3(Fx) — c^x. 
If a; is a feasible point, then jx is also feasible for all 7 S F. So the barycenter /3{Tx) is a 
convex combination of feasible points. The claim follows as the feasible region is convex. □ 

Since we assumed that LP(yl, b, c) has an optimal solution, the following is an immediate 
consequence of the preceding result. 

Corollary 14. There exists an optimal solution of ljP{A,b, c) which is a fixed point with 
respect to the entire automorphism group Aut(LP(v4, 6, c)). 

Up to minor technical details Theorem 3.3 of [S] generalizes Corollary [T3] to semi-definite 
programming. 

Let LP(A, 6, c) be a linear program with P(A,b) bounded and full-dimensional, and let 
r = (717 72, ■ • ■ , 7t) be a finite subgroup of Aut (LP {A, b,c)). Following Remark |9] we can 
compute a matrix E such that the kernel {x \ Ex = 0} is the fixed space Fixr(R"): for each 
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7i we determine a dual basis for the eigenspace {x \ (7^ — id)a; = 0} by solving a square 
system of linear equations. The total number of operations to do so is of order 0{tn^). 
Throughout this paper we measure algorithmic complexity in the RAM model; that is, we 
ignore the encoding lengths of real numbers, and all arithmetic operations are assumed to 
take constant time. The group F acts on the rows of the extended matrix (A|&), and we 
define a new extended matrix by summing the rows of the same F-orbit. We have 

the following general result. 

Theorem 15. The polyhedron 

P' = {xeW \ A'x < b', Ex = 0} 

is the set ¥ixr{P{A, b)) of feasible points which is fixed under the action o/F. In particular, 
P' = {P{Tx) I X G P{A, b)}. Each optimal solution of the linear program 

s.t. yE jx< yo j , X e 

is an optimal solution ofhP{A,b,c), and the objective values are the same. 

Proof. We constructed the matrix E to guarantee that each fixed point in P = P{A, b) 
satisfies the equation Ex — 0. Further, each inequality of the system A'x < 6' is a positive 
linear combination of valid inequalities. It follows that Fixr(P) is contained in P' . 

To prove the reverse inclusion consider a point x which is fixed by each transformation 
in F but which is not contained in P. Then for some index i we have the strict inequality 
ai^.x > bi. Without loss of generality we can assume that the first k rows ai,. , 02,., . . . , Ok,- 
of A form the F-orbit of the row ai^.. It follows that bi — b2 = ■ ■ ■ = bk = bi. Moreover, 
since x is a fixed point we have 

ai,.x = 02, .X = ■ ■ ■ — ai^.x = Oi^.x > bi . 

This implies that o,j.-)x > kbi, and hence x is not contained in P' . We conclude that 

P' is the set of points in P fixed by each transformation of F. Now Lemma [IT] says that 
P' is the image of P under the map x t-^ /3{Tx). The claim about the linear program (|3]) 
follows from Corollary [TH □ 

Remark 16. It has been observed by Scharlau and SchiirmanrQ that the vertices of the 
polyhedron P' are barycenters of orbits of vertices of P. This is a consequence of the fact 
that P' is the image of P under the linear map x n> /3(Fx). 

Corollary [M] and Theorem [15] yield a direct algorithm for solving a symmetric linear 
program: instead of solving LP (A, b, c) one can solve the linear program ((3|). The benefit is 
the following: The larger the group F < Aut(LP(A, &, c)) the smaller the dimension of the 
fixed space and the number of constraints. 

Remark 17. Formally, the feasible points of the derived linear program live in the same 
space M" as the original linear program. However, an algorithm based on the Simplex 
Method directly benefits if the solutions are contained in a proper subspace: the rows of 
the matrix E describing the fixed space never have to be exchanged in a Simplex tableau. 
Alternatively, one can project Fixr(R") onto a full-dimensional coordinate subspace, solve 
the projected linear program and lift back. 
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In the special case where the hnear program adniits a group of automorphisms acting on 
the standard basis of (that is, the groups acts by permuting the columns) it is standard 
optimization practice to identify variables in the same orbit, and to solve the reduced linear 
program. Theorem 1151 generalizes this approach to arbitrary groups of automorphisms. 

3. Symmetries of Integer Linear Programs 

We now turn to our main focus. Associated with LP(A, b, c) is the integer linear program 

, , max c*x 

^ ' s.t. Ax <b, xeZ" , 

which we denote as ILP(A, b, c). Throughout we make the same assumptions as above: the 
linear program LP(A, 6, c) is feasible, the matrix A does not have any zero rows, and the 
inequalities bijectively correspond to the hyperplane arrangement 6); see Remark [S] 

Definition 18. A symmetry of the integer linear program ILP(A, 6, c) is an automorphism 
of LP(y4, b, c) which acts on the signed standard basis {±ei, ±62, . . . , ±e„} of R" as a signed 
permutation. 

The symmetries of the integer linear program Q form a group Sym(ILP(A, b, c)) which 
is a subgroup of the group 0„Z, the group of all 0/1/— 1-niatrices with exactly one non-zero 
entry per row and column. We have 0„Z = 0„R fl GL„Z, and 0„Z is isomorphic to the 
Coxeter group of type i?„, the group of automorphisms of the regular ri-dimensional cube 
and its polar, the regular n-dimensional cross polytope. As a consequence, the group of 
symmetries of an integer linear program is finite, even if Aut(LP(yl, 6, c)) is infinite. 

The motivation for our definition is Lie-theoretical: let F be any finite subgroup of GL„Z. 
Then P is a compact subgroup of GL„M, hence it is contained in (a conjugate copy of) the 
maximal compact subgroup 0„R. It follows that, up to conjugation in GL„R, the group P 
is a subgroup of 0„Z. 

As an abstract group 0„Z is isomorphic to the wreath product 

Z2 I Sym(n) = (Z2)" >^ Sym(n) , 

where Z2 is the cyclic group of order two and Sym(n) is the symmetric group of degree n; 
the group Sym(n) acts on the direct product (Z2)" by permuting the factors. Each element 
of OnZ can be written as a product of a sign vector and a permutation. Since a permutation 
is a product of disjoint cycles, each signed permutation is a product of signed cycles which 
are disjoint. In terms of notation we write the signs between the indices within a cycle. This 
is to say, (1+2—4+3—) denotes the signed permutation matrix 



^0 





-1 




1 




















1 


^0 


-1 





0/ 



which is to be multiplied to column vectors from the left. 

Remark 19. In the optimization literature the authors often restrict their attention to 
symmetries permuting the standard basis vectors; for instance, see Margot and the 
references listed there. However, our more general analysis below shows that taking signed 
permutations into account does not cause any extra effort. Moreover, if the polyhedron 
P{A, b) is full-dimensional and bounded the group of automorphisms of the linear relaxation 
is already finite by LemmajT) Then Aut(LP(A, b, c)) n GL„Z is already contained in 0„Z by 
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the Lie-theoretical argument given above. Hence, at least in this case, considering groups 
of signed permutations is a natural choice. 

Before we will inspect groups of symmetries of integer linear programs we need to collect 
a few basic results on the action of the group 0„Z on the entire space M". Throughout let 
r be a subgroup of 0„Z. Then F acts on the standard basis 

S = {±ei, ±62, . . . , ±e„} . 

In the sequel we will always consider this particular action of F. There are two kinds of 
orbits to distinguish: the bipolar orbits contain at least one pair ±ei, while the unipolar 
orbits do not. Since F is a linear group, a signed permutation a € T with aci = eej and 
e € {±1} maps —Ci to —ecj. Hence, a bipolar orbit only consists of pairs, that is, — O = O. 
On the other hand, for each unipolar orbit O the set —O = {— e,; | € O} forms another 
orbit, and F acts equivalently on O and — O. 

Proposition 20. For the fixed space of F we have 

Fixr(M") = \m{l3{0) \0 orbit o/F} = lin{/3(0) \0 unipolar orbit o/F} . 

Proof. The first equality is a consequence of Lemma [T^ The second equality holds as 
(3{0) = for any bipolar orbit O. □ 

Remark 21. The points in S are the vertices of the regular n-dimensional cross polytope. 
If O C 5* is a unipolar F-orbit, then f3{0) is the barycenter of the non-trivial face of the 
cross polytope which is spanned by the vertices in O. In view of cone polarity the action of 
F on 5 is dual to the induced action on the vertices of the regular cube [—1,1]. That is, the 
two corresponding representations of F, on M" and its dual space, form a contra-gradient 
pair. 

We call the action of F on the set of signed standard basis vectors semi-transitive if there 
are precisely two opposite orbits of length n. Moreover, we call the action sub-transitive if 
there is no proper coordinate subspace which is invariant under F. Clearly, a semi-transitive 
action is necessarily sub-transitive. The converse does not hold, but we have the following 
characterization. 

Proposition 22. Suppose that F acts sub-transitively. Then either F acts semi-transitively 
with orbits O and —O such that the fixed space 

Fixr(M") = lin(/3o) = \m{P-o) 
is one-dimensional, or Fixr(K") — 0. 

Proof. If F has a bipolar orbit O, then O equals the entire set S of signed standard basis 
vectors because F acts sub-transitively. In this case the fixed space reduces to the origin. If, 
however, each orbit is unipolar, we have exactly one pair (O, —O) of opposite orbits, again 
due to sub-transitivity. Now the claim follows from Proposition [20l □ 

Corollary 23. IfT acts semi-transitively, then F is conjugate to a subgroup o/ Sym(n) in 
0„Z. 

Proof. Let O and — O be the two orbits of F, both of which have length n. Pick a transfor- 
mation e G 0„Z which maps O to the standard basis {ei, 62, . . . , e„}. Now for each 7 G F 
the conjugate transformation 676^^ leaves the sets {ei, 62, ... , e„} and {— ei, —62, . . . , — e„} 
invariant. We conclude that eFe""'^ is a subgroup of Sym(n). □ 
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We now interprete the results above for integer linear programming. Consider an integer 
linear program ILP(A, b, c) such that the set P{A, b) of feasible points of the linear relaxation 
is full-dimensional. Let T < Sym(ILP(A, &, c)) be a group of automorphisms. We have 
r < 0„Z. The action of F on the set {±ei, ±62, . . . , ±e„} can be decomposed into orbits. 
In this way the most relevant case occurs when F acts sub-transitively. From Lemma [7] we 
know that c is contained in the fixed space Fixr(ffi"), and then Proposition [21] says c 7^ 
enforces the action of F to be semi-transitive. Finally, by Corollary [23] we can conjugate F 
into a subgroup of Sym(n) acting on the standard basis {ei, 62, . . . , e„}. This is the situation 
that we will be dealing with in our algorithms below. 

4. Layers of Integer Points 

Our goal is to describe an algorithm for the efficient solution of a highly symmetric integer 
linear program. Again we consider ILP(A, b, c) with a group F of automorphisms as above. 

Let us assume that the objective function c ^ is projectively rational. This means 
that we require c to be a constant real multiple of a rational vector. For such a vector c let 
coprime(c) be the unique integral vector with coprime coefficients such that c = /3Coprime(c) 
for some positive real p. If c is a multiple of a standard basis vector, the single non-zero 
coefficient of coprime(c) is defined to be ±1. For an integer k the k-th c-layer is the affine 
hyperplane 

k 

Hc,k = kcr (.T M> c*x) H ^ coprime(c) . 

||coprime(c) II 

We have iJc,fc — Hpc^k for all p > 0, and H^c.k = ^Hc,k = H^^-k- All points in Hc^k attain 
the same value k with respect to the rescaled objective function coprime(c). We call k the 
number of the c-layer i?c,fc- The intersection of Hc^k with the line Rc is called the center. 

Lemma 24. If c ^ is projectively rational, the integral point x G Z" is contained in the 
c-layer with number coprime(c)*a;. 

Proof. The number k — coprime(c)*a; is an integer. We abbreviate d — coprime(c) and 
compute 

J. / k \ + / d xd \ d d J. J. 
cM Trd) = cU 7t) = ttCx = ex. 

Hence x — (k/ \\d\\^)d is contained in the kernel of the linear form c*, that is, the point x lies 
in the affine hyperplane H^.k- 

For the following result it is crucial that the coefficients of coprime(c) are coprime. 

Proposition 25. If c ^ is projectively rational, the c-layers Hck for A: G Z partition the 
set Z" of all integral points. 

Proof. From Lemma [U] is clear that each integral point is contained in some c-layer. By 
construction it is also obvious that the c-layers are pairwise disjoint. It remains to show 
that Hc^k n Z" is non-empty for all A; G Z. 

Let d = coprime(c) . Since the coefficients di,d2, . ■ . ,dn are coprime there are integral 
coefficients a;i, X2, . . . , x„ such that 

Xidi+X2d2 + V Xndn = gcd(di , ^2 , ■ • • , dji ) = 1- 

However, the left side of this equation equals c^x, whence the point x is contained in the 
first c-layer -ffc,i- Now c*{kx) — kc*{x) = k implies that the fc-th layer contains the integral 
point kx for arbitrary k & "L. □ 
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Input; {A, b) such that Syni(ILP(A, 6, 1)) acts transitively on standard basis 
Output: optimal solution of ILP(yl, b, 1) or "infeasible" 

1 let z = CI be a symmetric optimal solution of the LP relaxation LP(v4, b, 1) 

2 k ^ [n(\ 

3 repeat 

4 if it exists then 

5 I let a: be integral point in P{A, b) n Hi^k 

6 else 

7 \_ 

8 until feasible x found or A: < n[(\ 

9 return x or "infeasible" 

Algorithm A: Reduction to 1-layers 

Another way of putting the statement above is that coprime(c) is the unique generator 
of the unique minimal Hilbert basis of the one-dimensional pointed cone M>oc. 

Remark 26. An important consequence of Proposition [55] is that for any given bounds 
£,u £ there are only finitely many c-layers with feasible integral points whose objective 
values lie between i and u. This does not hold if the objective function is not projectively 
rational. 

Theorem 27. For given A and b such that Sym(ILP(A, b, 1)) acts transitively on the stan- 
dard basis the Algorithmic solves the integer linear program 1LP(A, 6, 1). 

Proof. Recall that throughout we assumed that the set of feasible points of the linear relax- 
ation is bounded. Hence it cannot occur that the integer linear program is unbounded. 

Let r < Sym(ILP(A, 6, 1)) be a transitive group of automorphisms. The fixed space is 
spanned by 1. If z is an optimal solution of the relaxation LP {A, 61), then, by Proposition[T3l 
the barycenter l3{Tz) = ^1 for ( = l/n(zi + Z2 + ■ ■ ■ + Zn) is also an optimal solution. 
Now [CJ 1 is an integral point in the fixed space with an objective value not exceeding the 
optimal value of the linear programming relaxation. Each 1-layer with a feasible integral 
point meets the one-dimensional polyhedron P' — {(3{rx) \ x £ P{A, b)}. We infer that no 
integral optimal solution of ILP(A, b, 1) can have an objective value strictly less than n[(\ . 

Due to Proposition [25] the 1-layers partition Z", and so the feasible points of ILP(A, 6, c) 
are contained in the set 

L"CJ 

U ^i.'^- 

fc=nLCJ 

□ 

The benefit of Algorithm [X] is that it reduces a (symmetric) n-dimensional integer linear 
programming problem to n integer feasibility problems in one dimension below. Since 
the latter is still an NP-complete problem not much is gained, in general. The situation 
changes, however, if we assume higher degrees of transitivity for the action of the group of 
automorphisms . 

Remark 28. Searching a family of parallel affine hyperplanes for integer points as in Algo- 
rithm [A] also plays a key role in Lenstra's algorithm for integer linear programming which 
requires polynomial time in fixed dimension [11) . 
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5. Searching Integer Layers Efficiently 



The question remaining is how to test ILP-fcasibihty of a c-layer in an efficient way. Our 
key observation is that some optimal integral solution is close to the fixed space if the group 
of symmetries acts sufficiently transitive. 

Definition 29. Given a c-layer with center z, an integral point in the c-layer is a core point 
if it minimizes the distance to z. 

Example 30. For the objective function c = 1 and an integer k = qn + r with g G Z and 
r G {0, 1, . . . , n — 1}, the set of core points in the k-th layer consists of all integer points with 
r coefficients equal to g'+l and n — r coefficients equal to q. In particular, the number of core 
points in this case equals (") . These core points are the vertices of an (r, n)-hypersimplex, 

translated by the vector ql. 

For the algorithms below the geometric structure of the set of core points is very relevant. 
We therefore make a short digression: the {r,n) -hyper simplex A(r, n) is the 0/1-polytope 
with vertices 



where 5* ranges over all r-element subsets of [n]. The hypersimplices are highly regular 

structures, and this yields the following. 

Proposition 31. Let T < GL„M be a linear group which acts jj -transitively on the standard 
basis. Then T acts transitively on the set of vertices of the {r,n) -hyper simplex for any 



Proof. By assumption F acts transitively on the r-element subsets of [n] for r < fi. Since F 
is a linear group it thus acts transitively on the set of vertices of A(r, n). The corresponding 
claim for the remaining hypersimplices follows since A(r, n) is affinely isomorphic to A(n — 



Below wo will apply the previous results in the special case where fj, > ln/2\. Then the 
groups acts transitively on the sets of vertices of all hypersimplices. 

Lemma 32. Let x G P{A, b) be an LP-feasible point in the k-th c-layer, and let 7 G 
Sym(ILP(A, 6, c)) with jx ^ x. Then any point in the interior of the line segment [x, 7a;] is 
LP-feasible and closer to the center of the k-th c-layer than x. 

Proof. Since 7 is an orthogonal linear map it preserves distances. The center z of the fc-th 
c-layer is fixed under 7, and this implies that {x, z,^x) is an isosceles triangle. We infer that 
\\p — z\\ < \\x — z\\ for all points p in the interior of [a;, 7a;]. Since 7 is an automorphism of 
the linear relaxation LP{A,b,c) the point 7a; is feasible, too. The feasible region is convex, 
and hence p is feasible. □ 

Theorem 33. Suppose that F < Sym(ILP(A, &, 1)) acts {[n/2\ -\- 1) -transitively on the 
standard basis o/M", and n>2. Then either each core point in the k-th 1-layer is feasible 

or Hi J, does not contain any feasible point. 

Proof. Let x he a feasible integer point in the k-Va 1-layer which is not a core point. We 
will show that there is another feasible integer point which is closer to the center, and this 
will prove the claim. 

Due to the invariance of Z" under translation by integer vectors we may assume that 
k G {0,l,...,n — 1}. Since x is not a core point, in particular, it is not the center of 




re {0,1,..., 



/i} U {n — /U, n — /X -|- 1, . . . , n} . 



r, n) via the map a; i->- 1 — a;. 



□ 
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the fc-th layer. Hence x is not contained in the fixed space Kl, which means that not all 
coordinates of x are the same. We split the set [n] of coordinate directions into two subsets 
by considering 

{i I Xi is even} and {i \ Xi is odd} . 

Then one of the sets — denoted by / — contains at least [(n + 1)/2J elements, while 
the other set J has at most [n/2j elements. We will employ the [n/2j -transitivity of the 
automorphism group to control J, and the additional degree of freedom to produce two 
distinct feasible integer points. We distinguish two cases. 

(i) Suppose that x has two different coordinates, say Xu and x^, which are in the 
same congruence class modulo two. That is, the set {m, w} is contained in either 
/ or J. Observe that this condition is satisfied whenever x has at least three 
pairwise distinct coordinates. Due to the {\n/2\ + l)-transitivity of F there is an 
automorphism 7 G F which leaves J invariant and which maps u to v. Since J 
is invariant, its complement I — [n] \ J is invariant, too. Notice that we do not 
require the set J to be non-empty (if {u,v} C /). 

Letting x' = 7a; we observe that Xi and x'^ are congruent modulo two for all 
i G [n]. Since Xu ^ Xv — x^(^u) = x'^ we have x ^ x' , and hence 

V = ^(a^ + a;') = ^(2^ + 72;) 

is an integer point in the interval [x,7x]. 

(ii) Otherwise the point x has exactly two different coordinates x^ and Xy , one of them 
being even, the other one odd. Without loss of generality, Xi = Xu for all i S / and 

for all j G J. Due to the transitivity of F there is an automorphism 7 G F 
with 7e„ — By. Then x and "fx are distinct points. Consider an interior point 

(5) y ^ Ax + (1 - X)jx for < A < 1 

in the line segment [x,7a;]. We want to find a parameter A such that y is integral. 
As X has only two distinct coordinates the i-th coordinate of y can attain the 
following values only: 

or 



or 



Since x is integral coordinates of types ([6]) and ([7|) are integers for arbitrary param- 
eters A G (0, 1). The coordinates of types ([8|) and ^ are integral if A - \xu —Xy \ G Z. 

We assumed that x is contained in the fc-th 1-layer for some fc = 0,l,...,n— 1 
and that it is not a core point. In Example [30] it has been observed that the core 
points in these layers are the vertices of a translated hypersimplex. We learned 
that some coordinate difference \xi — Xk \ must exceed one. Since all coefficients are 
equal to either Xu or Xy it follows that \xu — Xy \ > 2. We can now set 



(6) 




— A.T„ - 


Ml 


- X)Xy 




(7) 


Vi 


— Xxy - 


Ml 


- X)Xy 


— IC^ 


(8) 




— XXy - 


Ml 


— X)Xy 


= A(x- 


(9) 




= XXy - 


Ml 


- X)Xu 


= X{x- 



X = -± 



in the formula ([5]). 
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Input; {A, b) such that Syni(ILP(A, b, 1)) acts ([n/2j + l)-transitively on standard 
basis 

Output: optimal solution of ILP(A, 6, 1) or "infeasible" 

1 let z = CI be a symmetric optimal solution of the LP relaxation LP(74, b, 1) 

2 d [nCJ - n [CJ 

3 repeat 



4 



X 



(LCJ+i,...,LCJ+i,LCJ,.-.,LCJ) 



d n — d 

if X infeasible then 
\_d^d-l 

7 until feasible x found or d < 

8 return x or "infeasible" 



Algorithm B: Core point algorithm 



In both cases we obtain an integral point y in the interior of the interval [x,7x]. By 
Lemma [32l such a point is always closer to the center than x. This shows that there exists 
a feasible core point in the same layer as x. Applying Proposition [311 with ji = [n/2\ + 1 
yields that then each core point must be feasible. □ 

Now Algorithm |X] can be modified in Step [5] to check a single core point per layer for 
feasibility, provided that the group of automorphisms of the ILP acts at least ([?t./2J + 1)- 
transitively. This is our Core Point Algorithm [Bl 

Corollary 34. For given A and b such that Sym(ILP(yl. b, 1)) acts {[n/2\ + 1) -transitively 
on the standard basis the Core Point Algorithm\B\solves the integer linear program ILP(A, b, 1) 
in 0{mn^) time. 

Proof. The correctness follows from Theorems [57] and The main loop of the algorithm 
is executed at most n times. In each step the costs are dominated by checking one point in 
for feasibility against m linear inequalities. □ 

Remark 35. The linear search in Algorithms 1X1 and [B] cannot be substituted by a direct 
bisectional approach. The reason is that the set of all fc in {0, 1, ... , [n^J — n[CJ} such that 
the fc-th 1-layer contains a feasible point is not necessarily (the set of integer points of) an 
interval. 

6. Finding All Symmetries 

For the algorithms presented it is never necessary to know the entire group of automorphisms 
of LP{A,b,c) or ILP(A, 6, c). Generally, any subgroup will do, the larger the better. Yet 
here we would like to discuss the question of how to find automorphisms of integer linear 
programs. From the input data we will construct a labeled graph G{A, b, c) whose group of 
labeled automorphisms coincides with Sym(ILP(yl, 6, c)). 

Expressing symmetry in optimization via graph automorphisms is not a new idea: the 
linear automorphism group of a polytope and of a linear program can be obtained by com- 
puting the automorphism group of a certain graph as described by Bremner, Dutour Sikiric, 
and Schiirmann ;4j. The combinatorial automorphisms of a polytope are the (labeled) graph 
automorphisms of the bipartite graph encoded by the vertex-edge-incidences. This directly 
follows from the fact that the face lattice of a polytope is atomic and coatomic; see Kaibel 
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and Schwartz [S] . Liberti studies automorphisms of optimization problems which are more 
general than integer linear programs |12| . His approach, however, deals with expression 
trees obtained from a specific encoding of the optimization problem. None of these concepts 
seems to be directly related to the kind of symmetry studied here. An idea similar to ours, 
however, has been applied by Berthold and Pfetsch [T] to find symmetries of 0/1-ILPs. 

The complexity status of the graph isomorphism problem is notoriously open. While the 
known algorithms for determining the automorphism group of a graph require exponential 
time, there exist software packages, for instance, nauty |16j or SymPol [211, that can solve 
this problem very well in practice. 

For a given matrix A G M™^", right hand side b G M™, and objective function c G R" we 
will now associate two undirected simple graphs, the ILP graph G(A, 6, c), and the restricted 
ILP graph G'{A, b, c). For the sake of a simplified exposition we start out by describing the 
restricted ILP graph. Throughout we assume that the rows of the extended matrix are 
normalized as described in Remark [5] We have one node for each position in the matrix 
A, one node pi for each row, and one node Cj for each column, that is, G [m] x [n], 

where [n] = {1, 2, . . . Further, we have one node k„ for each distinct coefficient u in 
the matrix A, one node A,, for each distinct coefficient v of 6, and one node for each 
distinct coefficient w of c. This gives a total of mn + m + n + ua + rib + ric nodes, where 
riA, rib, and ric denotes the respective number of different entries in A, b, and c. The nodes 
receive labels in the following way: all positions share the same label, the rows receive a 
second, and the columns a third label. Each node corresponding to one of the coefficients 
receives an individual label. This way we arrive at ua + rit + ric + 3 labels altogether. The 
edges of G'{A,b,c) are defined as follows: the node a^- is adjacent to pi and (j as well as 
to the coefficient node which represents the coefficient aij of the matrix A. Moreover, the 
row node pi is adjacent to the node A^. , and the node Q is adjacent to the node fic^. This 
totals to 3mn + m + n edges. 

Example 36. The reduced ILP graph of the integer linear program 

< 3 

< 3 

< 3 , Xi gTL 
is shown in Figure [TJ 

Let 7 be an automorphism of G'(v4, 6, c) which respects all node labels. Since the common 
label of the column nodes is preserved 7 induces a column permutation i]}^ S Sym(n). Now 
'07 acts on the standard basis {ei, 62, . . . , e„}, and by linear extension we obtain a linear 
transformation which we denote V'^ • 

Lemma 37. The linear transformation ip* is a symmetry 0/ ILP(^, 6, c). 

Proof. As above let 7 be a labeled automorphism of G' = G'{A, b, c) with induced column 
permutation ip — tp-y and linear transformation ip* € Sym(rt) < GL„R. As for the column 
nodes the graph automorphism 7 also induces a permutation (f) S Sym(m) of the row nodes 
of G'. The position nodes a-ij form a label class of their own, and so they are permuted by 
7 as well. Since each position node is adjacent to precisely one row and one column node 
we infer that 7(0;^) = CK0(i),i/i(j) ■ Each position node is adjacent to precisely one matrix 
coefficient node, each of which forms a singleton label class. This implies that the coefficient 
aij corresponding to the node aij is the same as the coefficient a^i^i-^ ipi^jy Likewise we obtain 
bi = b^(i) and Cj = c^(j) ■ This means that ip* is a symmetry of ILP(A, 6, c). □ 



max xi + X2 + X3 

, - s.t. xi + 2x2 

^ ^ X2 + 2x3 

2xi + 
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Proposition 38. The map j ^ tp* is an isomorphism from the group of labeled automor- 
phisms of the graph G' {A, 6, c) to the group Sym(ILP(A, b, c)) H Sym(n). 

Proof. We describe the inverse map. To this end let a he a symmetry of ILP(^, 6, c) which 
acts on the standard basis of M". Hence a induces a permutation (j> of the rows of the 
extended matrix {A\b) and a permutation ^p of the columns of A. It is obvious how 4> and 
"0 induce permutations of the row nodes and of the column nodes of G'. By the same 
reasoning as in the proof of Lemma l37l the pair (0, V') uniquely extends to a labeled graph 
automorphism 7(ct) of the reduced ILP graph. 

We omit the straightforward proofs that the equations ^{4^^) — 7 and V^^^^) ~ both 
hold. From these it follows that the map 7 H> is bijective. In both groups the multipli- 
cations are given by concatenations of maps. A direct computation yields ~ ''/''7i^*2' 
all maps are acting on the left. Hence the group structures are preserved. □ 

We now explain how the full ILP graph G{A, b, c) differs from the restricted ILP graph 
G'{A,b,c). The key to the construction of G'{A,b,c) is the map 7 i-^ V'^ yielding a linear 
transformation which acts as a permutation of the standard basis of R". In order to allow for 
signed permutations certain nodes have to be duplicated: each column node Q in G'{A, 6, c) 
gets a twin node Q in G{A,b,c), each matrix coefficient node corresponding to a non- 
zero coefficient gets a twin node aij. Moreover, we add further nodes representing negatives 
of non-zero coefficients in the matrix A and the objective function c unless nodes with 
these labels already exist. This way ILP (A, 6, c) has less than twice as many nodes as 
ILP' (A, b, c); it is always strictly less as the nodes corresponding to the coefficients in b are 
never duplicated. We also add edges such that first atj is adjacent to pi and for all i 
and J, second Cj is adjacent to /^-c^, third aij is adjacent to K-a^j, and, finally, the twins 
are matched up: aij is adjacent to ctij and Q is adjacent to Q. The labeling is extended 
in a way such that twins share the same label; the nodes newly introduced for negatives of 
coefficients receive new singleton labels. 

Each labeled graph automorphism of G'{A,b,c) uniquely extends to a labeled graph 
automorphism of G{A, b, c), but the automorphism group of the non-reduced ILP graph is 
larger, in general. We have the following result. 
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Theorem 39. The group of labeled graph automorphisms of G{A,b, c) is isomorphic to the 
group of symmetries ofILP{A,b,c). 

Proof. One can follow the strategy in the proof of Proposition [38l We know that a labeled 
graph automorphism of G'{A, b, c) encodes a symmetry of ILP(yl, b, c) which permutes the 
set {ei, 62, ... , e„}. Now a labeled graph automorphism of G{A, 6, c) may map a column 
node Cj" to some node C^.. But then it follows that Q is mapped to since Q is the only 
column node adjacent to and Cfc is the only column node adjacent to Ck- This shows that 
the permutation of the column nodes can be extended to a linear transformation. As in the 
proof of Proposition [38] one can show that this linear transformation is a symmetry of the 
integer linear program. Conversely, each such symmetry acts like a signed permutation on 
the signed standard basis and yields a labeled isomorphism of the graph G(^, b, c). □ 

Roughly speaking, a class C of graphs is graph isomorphism complete if the problem of 
deciding isomorphy for any two graphs in C is as difficult as for general graphs, up to a poly- 
nomial time transformation. For a precise definition, for instance, see the monograph (lOj . 
The next result is not only of theoretical interest. To the contrary, for practical applica- 
tions it can be read as: finding the symmetries of an integer linear program via reducing to 
automorphisms of suitable (labeled) graphs, is the right thing to do. 

Theorem 40. The classes of ILP graphs and reduced ILP graphs are both graph isomor- 
phism complete. 

Proof. We only prove that the class of reduced ILP graphs is graph isomorphism complete. 
It is known that the class of bipartite graphs is graph isomorphism complete. Hence it 
suffices to encode an arbitrary bipartite graph as a reduced ILP graph, which is not too 
large. 

Let G ~ {V, E) be an undirected bipartite graph with m + n nodes V = U U W = 
{ui, . . . ,Um} U {wi, . . . ,Wn}. As our matrix Ac = (aij) G M™^" we take the bipartite 
adjacency matrix of G, that is. 



1 otherwise . 

For a second bipartite graph G" it is easy to see that the reduced ILP graph of ILP(j4g, 1, 1) 
is isomorphic to the reduced ILP graph of ILP(^g/, 1, 1) if and only if G is isomorphic to 



Remark 41. Rehn investigates arbitrary automorphisms of the integer lattice Z" in the 
context of polyhedral geometry [20] . In particular, his modification of a backtracking al- 
gorithm of Plesken and Souvignier |18j allows to obtain matrix generators of the group of 
symmetries. For practical applications this should be superior to our approach via graph 
automorphisms if the number m of constraints is much larger than the dimension n. 



In this section we will construct a specific class of highly symmetric convex polytopes among 
which one can find examples of rather high Gomory-Chvatal rank. The motivation for this 
construction is rooted in the systematic study of Gomory-Chvatal cuts and cutting-plane 
proof systems. Pokutta and Stauffer |19| propose a new method for computing lower bounds 
on the Gomory-Chvatal rank for polytopes contained in the 0/1-cube, and the polytopes 
constructed here provide examples which asymptotically almost attain the bounds obtained. 




G'. 



□ 



7. Hypertruncated Cubes 
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The subsequent section on computational experiments also contains results about these 
polytopes. 

Our construction starts out with the unit cube C — [0, 1]". Intersecting C with the 
hyperplane defined by ^ = r for r S {2, 3, . . . , n — 1} gives the hypersimplex A(r, n) 
which already appeared in Example 1301 Here we are interested in the {r^n) -truncated cube 
C — {x e [0, 1]" \ J2^i — We make one more modification to the cube C by defining 
the polytope 

(11) C" = conv(C' U Al) for A > r/n 

which we call the {r,n; X)-hypertruncated cube. Notice that the full group Sym(n) acts on 
the cube C as well as on the truncated cube C" as well as on the hypertruncated cube C". 
Hence our algorithms above can be applied. Our next goal is to describe the vertices and 
the facets of C". 

Proposition 42. Let n > 2, r G {2, 3, . . . , n — 1}, and r/n < A < 1. The vertices of the 
{r,n; X) -hypertruncated cube C" are 

es for all S C [n] with < r and Al . 

Proof. The points 65, for S C [n] and < r, are the vertices of the (r, n)-truncated cube 
C". They are also vertices of C" . Since nX exceeds r, the hyperplane = riA does not 

separate C" , and its intersection with C" only contains the point Al. Hence the latter is 
a vertex, too. Looking at the defining Equation (fTTj) shows that there cannot be any other 
vertices. □ 

Of course, the vertices determine the facets completely. In this case, it is particularly 
easy to read off the facets of C" by looking at the facets of C" and analyzing what changes 
in case the point Al is added as a generator. This proves the claim in [IQl Remark 3.3]. 

Proposition 43. Let n >2, re {2, 3, . . . , n — 1}, and X > r/n. The facets of the (r, n; X)- 

hypertruncated cube C" are 

(12) X, > , X, < I 

(13) (l-n+^)a;, + Vxfc > r 

(14) (l-r + A(n-l)).T, + (1-A)^.T;,. < X{n ~ r) 
for i G [n]. Ln particular, C" has precisely An facets. 

Proof. The facets of type (fT2| are the facets of the unit cube C. Together with the truncating 
inequality '^Xi < r they also form the facets of the truncated cube C . The remaining facets 
of C" are the facets through the vertex Al. Each of them is the convex hull of Al and a 
ridge of C contained in the truncating facet. A ridge is a face of codimension 2, that is, a 
facet of a facet. As pointed out above the truncating facet is the hypersimplex A(r, n). Its 
facets arise from the intersection with the cube facets. A hypersimplex facet of type {x G 
A(r, n)\xi = 0} is a deletion facet, and a hypersimplex facet of type {x G A(r, n) \ Xi = 1} 
is a contraction facet. The n — 1 points re^ for k ^ i span an in — 2)-dimensional affine 
subspace A containing the i-th deletion facet. However, these points are not contained in 
A(r, n). Looking for an affine hyperplane containing A and Al results in a rank-1 system 
of linear equations. This way we obtain the n linear inequalities of type (fT3|) . Similarly, the 
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afRne span of a contraction ridge is generated by the n — 1 points + (r — l)efe for k ^ i. 
Via the same approach we arrive at the n hnear inequahties of type (1141) . □ 

Remark 44. Pokutta and Stauffer [TH] show that the Gomory-Chvatal ranks of the (r, n; A)- 
hypertruncated cubes for r = [«/eJ, where e = 2.7172... is Euler's constant, and A = 
(m — 1)/to approach n/e — o(l) as to € N goes to infinity. In our experiments below we look 
at the case r — [n/ej and A = 1/2, that is, to = 2. 

8. Computational Results 

The following experiments were carried out on an Intel(R) Core(TM) 17 920 2.67GHz ma- 
chine, with 12GB of main memory, running Ubuntu 10.04 (Lucid Lynx). The performance 
of each core is estimated at 5346.16 bogomips each. All our tests were run single-threaded. 

The goal of the experiments is to compare the performances of a conventional branch- 
and-cut approach (with automated symmetry detection) and the Core-Point- AlgorithmlBlon 
highly symmetric integer linear programs. For the test of the conventional branch-and-cut 
method we used CPLEX, Version 12.1.0, while the Core-Point- Algorithm was implemented 
and tested in polymake, Version 2.9.9 [7]. As a major difference polymake employs exact 
rational arithmetic (via GMP), while CPLEX uses floating-point arithmetic. It should be 
stressed that CPLEX can detect if the symmetry group of an integer linear program contains 
the full symmetric group acting on the standard basis of R", and this is exploited in its 
algorithms. For input in this category (which includes all our examples below), it is thus 
quite a challenge to beat CPLEX. 

8.1. Hypertruncated Cubes. We tested our algorithms on the ( [n/eJ . n; l/2)-hyper- 
truncated cubes; see Remark 1441 In this case we have only 4n linear inequalities from Propo- 
sition|43]as input. Each coefficient is small, and computationally accuracy (for floating-point 
computations) or coefficient growth (for exact arithmetic) is not an issue here. This benign 
input can be dealt with easily up to high dimensions. Table [T] lists the timings for CPLEX' 
Branch-and-Cut and polymake's Core Point Algorithm. The timings required to obtain the 
solution of the linear relaxation are given separately for both systems. 

The fact that polymake takes more time is due to the overhead induced by the GMP exact 
rational arithmetic. Since coefficient growth does not occur the overhead versus floating- 
point arithmetic can be estimated to be constant. Hence the roughly quadratic overhead (in 
dependence of n) versus the CPLEX result is a consequence of the total algorithmic complexity 
of Oimv? ) from Corollary [Ml Altogether both solvers behave pretty well for these kinds of 
examples. 

An industry strength solver as CPLEX comes with a number of bolts and whistles which 
allow to tune its behavior in many ways. For the hypertruncated cubes this does not play 
any role. Since no parallel implementation of the Core Point Algorithm is available (yet) 
we set the number of CPLEX' parallel threads to one. 

8.2. Wild Input. One way to produce symmetric input for (integer) linear optimization 
algorithms is by brute force: One can take any system Ax < & of linear inequalities and let 
the full group Sym(?i) act. This way each original inequality may give up to n\ inequalities 
in the resulting symmetrized system. The symmetrized system of inequalities is Sym(n)- 
invariant by construction. In order to produce input to our algorithms which is less well 
behaved than the hypertruncated cubes studied above we will apply this procedure to a 
special class of polytopes, which can be considered "wild" . We aim at symmetric polytopes 
with many facets whose coordinates are not so nice, but still somewhat under control. 
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CPLEX polymake 



d 


time LP (s) 


time IP (s) 


time LP (s) 


time IP (s) 


100 


0.00 


0.07 


0.01 


0.08 


200 


0.07 


0.29 


0.02 


0.57 


300 


0.19 


0.73 


0.03 


1.88 


400 


0.41 


1.58 


0.06 


4.26 


500 


0.90 


2.99 


0.10 


8.39 


600 


1.52 


4.80 


0.14 


14.30 


700 


2.21 


7.01 


0.18 


22.29 


800 


3.44 


11.59 


0.24 


32.96 


900 


5.17 


16.37 


0.31 


47.11 


1000 


6.77 


21.66 


0.38 


65.11 


1100 


8.25 


26.55 


0.47 


85.52 


1200 


11.75 


35.47 


0.55 


111.59 


1300 


14.05 


45.63 


0.65 


142.76 


1400 


18.96 


57.23 


0.74 


175.23 


1500 


23.42 


73.50 


0.85 


217.19 


1600 


28.11 


78.13 


0.98 


263.21 


1700 


32.07 


97.23 


1.11 


315.38 


1800 


41.82 


128.88 


1.25 


374.63 


1900 


44.68 


137.22 


1.40 


444.48 


2000 


50.39 


154.35 


1.54 


511.59 



The first building block of our construction is the regular hexagon H whose vertices are 
at distance 56/6 from the origin, that is, 

f 56 

H = conv I —e'^"/^ 



A; = 0,1, 

D 

Notice that only in the formula above the letter 'i' denotes the imaginary unit, and we 
identify the complex numbers with M^. The coordinates of H are irrational; however, the 
subsequent steps in the construction are chosen such that we will arrive at a rational polytope 
in the end. The second item is the regular cross polytope scaled by 73/10, that is, 

f 73 

C{d) = conv < ^Yq^^ i & [d] 

Finally, we consider the join P * Q oi two polytopes P C and Q CW-, which is defined 
as 

P*Q = conv( {(a;, 0, 1) G X X M \ xeP} 

u{(0,?/,-l) e M'' X X R \ yeQ}). 

If P and Q arc full-dimensional polytopes with ji and v vertices, respectively, the join P *Q 
has dimension (5 + e + 1 and ji + v vertices. For the combinatorics of P * Q the exact values 
for the (5 + e + l)st coordinate are inessential, as long as they are distinct. We now replace 
the "—1" for the second factor by —11/12 to obtain the distorted join 

J{d) = conv( {(a;, 0, 1) G X M'* X M \ x & H) 

U j(0,y,-^) G X M*^ X M y&C{d) 



12' 
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Table 2. Symmetrized distorted joins of a hexagon with cross-polytopes 



CPLEX polymake 



d 


time LP (s) 


time IP (s) 


time LP (s) 


time IP (s) 


3 


0.00 


0.01 


0.00 


0.00 


4 


0.00 


0.06 


0.01 


0.00 


5 


0.00 


0.17 


0.01 


0.02 


6 


0.05 


0.74 


0.04 


0.04 


7 


0.13 


2.71 


0.09 


0.13 


8 


0.62 


10.15 


0.24 


0.38 


9 


2.08 


42.06 


0.69 


1.03 


10 


8.02 


135.51 


1.86 


2.89 



of the hexagon H with the cross polytope C{d). This polytope is further modified in two 
steps: First we perturb by rounding the (rational and irrational) coordinates to three decimal 
places and treating these as exact rational numbers. Since the polytopes H, C{d), and J{d) 
are simplicial this perturbation does not change the combinatorial types. Secondly, we 
symmetrize the polytope by letting the group Sym((i + 3) act on the facets of the perturbed 
polytope. The resulting inequalities form the input of our second class of experiments. 

The parameters 56/6, 73/10, and 11/12 which occur in the construction are chosen, 
more or less, at random. They do not have a specific meaning. We refrain from further 
investigating these symmetrized distorted joins and the geometry of lattice points inside. 
This would be tedious and at the same time irrelevant for our purposes. 

The interesting fact is that we get symmetric polytopes which are somewhat complicated, 
because they have lots of inequalities: for instance, yielding 885,768 inequalities for d = 10. 
As a consequence CPLEX cannot deal with these examples in a fully automated way. The 
best parameter settings that we found were 

parallel thread count: 1 

presolve indicator: no 

feas. pump heuristic: -1 

RINS heuristic: -1 

MIP optimization emph.: 2 

But even with these adjustments our implementation outperforms CPLEX by a large margin; 
see Table [21 This holds in spite of the fact that polymake computes with exact rational 
numbers throughout. 
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